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The ^^0(p,7)^'^F reaction rate is revisited with special emphasis on the stellar temperature range 
of T=60-100 MK important for hot bottom burning in asymptotic giant branch (AGB) stars. We 
evaluate existing cross section data that were obtained since 1958 and, if appropriate, correct pub- 
lished data for systematic errors that were not noticed previously, including the effects of coincidence 
summing and updated effective stopping powers. The data are interpreted by using two different 
models of nuclear reactions, that is, a potential model and R-matrix theory. A new astrophysical 
S-factor and recommended thermonuclear reaction rates are presented. As a result of our work, the 
^''0(p,7)^'^F reaction has now the most precisely known rate involving any target nucleus in the 
mass A > 12 range, with reaction rate errors of about 7% over the entire temperature region of 
astrophysical interest (T=0.01-2.5 GK). The impact of the present improved reaction rate with its 
significantly reduced uncertainties on the hot bottom burning in AGB stars is discussed. In contrast 
to earlier results we find now that there is not clear evidence to date for any stellar grain origin 
from massive AGB stars. 

PACS numbers: 



I. INTRODUCTION 

The "'^^0(p,7)"'^^F reaction is characterized by a number 
of exceptional attributes. At lower bombarding energies 
it provides a textbook example for a nonresonant reaction 
cross section since the lowest-lying resonance is located 
at a relatively high laboratory energy of 2.66 MeV jj. 
The absence of low-energy resonances and the high bind- 
ing energy of ^^O are among the main reasons for the fact 
that the cross section can be described in terms of simple 
nuclear reaction models. Indeed, the ^^0(p,7)^^F reac- 
tion is a prime example for the direct capture reaction 
model which assumes that the projectile is captured via a 
single-step process into a final state orbit outside a closed 
^^O core (see, for example, Ref. [I])- The absence of 
low-energy resonances is also the reason for the fact that 
the ^^0(p,7)"'^'^F reaction is the slowest process among 
all the proton-induced reactions in the CNO target mass 
region [3|. This reaction has been studied many times 
at low energies. The experimental techniques applied in- 
clude the activation method [1, [^, the in-beam study 
of prompt 7_-rays 0, 0, 13 and measurements in inverse 
kinematics |8|. It is generally assumed that the different 
measurements are in agreement. Despite these facts, the 
thermonuclear reaction rates evaluated by the NACRE 
collaboration have relatively large errors amounting, for 
example, to -1-35% and -43% at T=0.06-0.1 GK g. 

Many stars, including the Sun, will eventually pass 
through an evolutionary phase which is referred to as 



the asymptotic giant branch (AGB; see Ref. for a re- 
view) . This phase involves a hydrogen and a helium shell 
that burn alternately surrounding an inactive stellar core. 
Convection carries the products of nucleosynthesis to the 
stellar surface where material is ejected via strong stellar 
winds A fraction of the ejected matter condenses in 



form of small grains. Some of the grains travelled the in- 
terstellar medium and have been incorporated into prim- 
itive meteorites at the birth of the solar system. Such 
stellar grains can be extracted in the laboratory from 
their host meteorites and their isotopic composition can 
be analyzed with high precision [12] . The measured iso- 
topic abundance ratios deviate substantially from solar 
system abundance ratios, which result from a homoge- 
neous mixture of contributions from countless stars prior 
to solar system formation. Stellar grains, instead, exhibit 
isotopic abundance ratios that are characteristic for the 
composition of their parent stars, and hence they most 
likely provide strong constraints for stellar models. The 
present work has been motivated by studies showing that 
variations of the ^^0(p,7)^^F reaction rate influence sen- 
sitively the -'^^O/^^O isotopic ratio predicted by models 
of massive (> 4:Mq) AGB stars, where proton captures 
occur at the base of the convective envelope (hot bottom 
burning, HBB). Specifically, a recent study demonstrated 
that a fine-tuning of the ^^0(p,7)^^F reaction rate may 
account for the measured anomalous ^'^O/^^O abundance 
ratio in the extraordinary presolar spinel grain 0C2 p^ . 

In the present work we focus our attention on a de- 
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FIG. 1: Nuclear energy level diagram for ^^0(p,7)^'^F. The 
level parameters and the reaction Q-value are adopted from 
Ref. [3 and Ref. [3, respectively. 



tailed evaluation of existing ^^0(p,7)^^F cross section 
data (Sec. II). This level of detail was not practical 
in the broad NACRE compilation which includes many 
different reactions. We interpret the evaluated data in 
terms of two different nuclear reaction models (Sec. Ill) 
and extract astrophysical S-factors, reaction rates and 
associated errors (Sec. IV). The new ^^0(p,7)^^F ther- 
monuclear reaction rates are then incorporated into stel- 
lar models in order to derive improved oxygen isotopic 
ratios from the hot bottom burning in AGB stars (Sec. 
V). 



II. EVALUATION OF EXISTING DATA 

The radiative proton capture on ^^Q (Q=600.27±0.25 
keV [l3l) can proceed via El radiation either to the 
ground state (J'^=5/2+) or to the first excited state 
(E^=495.33±0.10 keV; J'^=l/2+) of ^^F A nuclear 
level diagram is shown in Fig. [T] In this section, we evalu- 
ate the existing cross section data and discuss certain cor- 
rections that we deem necessary. Notice that we make a 
distinction between experiments that measure separately 
the transitions to the two final states and those that re- 
port only on the sum contribution of both transitions. 
The data will be displayed in terms of the astrophysical 
S-factor, defined as 

S{E) = E a{E) e^'^" (1) 

with a the reaction cross section, e~^'^'' the Gamow fac- 
tor, and rj the Sommerfeld parameter. Throughout this 
experimental section, all kinematic quantities are given 
in the laboratory system, unless mentioned otherwise. 



A. Data of Chow, Griffiths and Hall (1975) 

The work of Chow, Grifhths and Hall reports on 
cross sections for ^^0(p,7)^^F at seven bombarding ener- 
gies in the range of Ep=845-2554 keV. The cross sections 
are separately given in their Tab. 4 for the transitions to 
the ground state (4 data points) and first excited state 
(7 data points) of ^^F. It is important to emphasize that 
these authors normalized the radiative capture cross sec- 
tion to the Rutherford scattering cross section at Ej,=405 
keV. This procedure provides capture cross sections inde- 
pendently from the target stoichiometry, stopping powers 
and beam straggling and thus avoids many normalization 
errors that are frequently associated with absolute cross 
section determinations. For more information, see Ref. 

i. 

Not all the experimental details are provided in Ref. 
Q in order to assess the quoted errors. Some of their re- 
ported cross section errors are as small as ±3%. However, 
the elastic scattering measurement alone contributes an 
error of about ±2%, as can be seen from their Tab. 1. 
Based on the geometry of the experiment and the meth- 
ods applied, the following estimates for additional error 
contributions seem reasonable: ±3% for the 7-ray effi- 
ciency; ±1% for the escape peak detection; ±1% for the 
angular distribution correction; and ±3% for the deter- 
mination of the effective bombarding energy in their 25 
keV thick solid oxygen targets. Other sources of errors, 
such as corrections for the finite solid angle of the 7-ray 
detector, are mentioned in Ref. ^61]. The quadratic addi- 
tion of the above mentioned errors gives a value of ±5% 
and, therefore, we suspect that the smallest of the cap- 
ture cross section errors quoted in their Tab. 4 have been 
somewhat underestimated. As a result of the above dis- 
cussion, we adopt errors of at least ±5% for their quoted 
capture cross sections. The modified data of Chow, Grif- 
fiths and Hall @ are displayed in Fig. 2. 

B. Data of Morlock et al. (1997) 

The measurement of Morlock et al. [tI (see also Ref. 
[Tst ) was performed in the energy range of Ep=280-3740 
keV by using an extended gas target. Their ^^0(p,7)^^F 
cross section was determined relative to the ^^0(p,p)^^0 
elastic scattering cross section which, in turn, was nor- 
malized to the Rutherford scattering of protons on heavy 
noble gases (krypton and xenon). A number of effects 
were considered by Refs. 0, HBl in the analysis, includ- 
ing the energy loss and straggling of protons in the gas 
target, the energy spread of the beam, the finite solid 
angle of the 7-ray detector, the angular distribution of 
the emitted 7-rays, and the dependence of the detection 
efficiency on 7-ray energy and emission angle. 

We noticed, however, that the data of Refs. 0, 
were not corrected for coincidence summing. This effect 
arises since the primary 7-ray populating the Ea;=495 
keV level and the secondary 7-ray from the subsequent 
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de-excitation of this state are in time-coincidence (Fig. 
1). Since both photons are emitted nearly simultane- 
ously, the detector is not able to distinguish between the 
two, but instead records a single event with the summed 
energy. This experimental artifact gives rise to too many 
counts for the ground state transition ( "summing- in" ) , 
and too few counts for the transition to the first excited 
^^F state ( "suniming-out" ) . The effect can only be cor- 
rected for with the proper knowledge of both the peak 
and total detection efficiencies (see Ref. [1] for details). 
Most of these quantities have not been reported in Refs. 
0) EBl, but can be deduced from other sources, as de- 
scribed below. 

We adopted the peak efficiencies from Ref. [l^ and 
from Refs. [lB)ll3|- The latter studies were performed 
with a similar detection setup as in Morlock et al. 
In order to estimate the total detection efhciency, we 
performed GEANT simulations for a similar setup in- 
volving a Germanium 7-ray detector of 100% relative 
efficiency. We find summing-out correction factors of 
1.06±0.02 for the primary transition to the first excited 
"'^'^F state and 1.05±0.02 for the secondary de-excitation 
of the 495 keV level. For the primary ground state tran- 
sition, the summing-in depends on the cross section of 
the transition DC^495. The correction reduces the re- 
ported cross section at the lowest measured energies (few 
100 keV) by about 10%, while the correction is negligible 
at the highest energies (several MeV). 

We found evidence in Morlock's thesis [l^ that an ad- 
ditional error of 10% has to be applied to all data points. 
Note that this additional error was disregarded in Ref. 
Q • We also discovered an inaccuracy (too few significant 
digits) in the numerical conversion from cross section to 
S-factor presented in Ref. [lB|- The modified data of 
Morlock et al. (3] and Morlock [TB| are displayed in Fig. 
2. 



C. Data of Hester, Pixley and Lamb (1958) 



with Zt the charge of the target nucleus, E the labora- 
tory bombarding energy in MeV, Y the thick target yield 
and ee// the effective stopping power. This expression is 
derived under the following assumptions: (i) the target 
is infinitely thick; (ii) the S-factor is approximately con- 
stant in the energy range of the measurement; (iii) the 
stopping power is constant over the effective thickness 
of the target. For the effective stopping power of their 
AI2O3 target, eg// = eo -l- (2/3) e^;, they used a value of 
Ee// = 2.8 X 10~^^ eVcm^ that was adopted from an un- 
published report l^l- We corrected the measured cross 
sections of Ref. ^ | by using modern stopping powers de- 
rived from the code SRIM [19,] . The correction decreases 
the cross section reported in Ref. [1] by 8% (2%) at the 
highest (lowest) energy measured. The modified data of 
Hester, Pixley and Lamb [3| are displayed in Fig. 2. 



D. Data of Becker et aL (1982) 

Becker et al. [sj measured the cross section for ra- 
diative capture on ^^O at a center-of-mass bombarding 
energy of Ep'"=853 keV. Their study was performed in 
inverse kinematics by directing a -'^^O heavy- ion beam 
onto an extended windowless hydrogen gas target. It 
must be emphasized that Becker et al. performed simi- 
lar measurements for other nuclear reactions (see Tab. 1 
in Ref. 0]) and that their reported resonance strengths 
are frequently used in the literature as absolute strength 
standards (see, for example, Refs. [1, [13, HH)- For the 
^^0(p,7)^''F cross section at Ep'"=853 keV they report a 
value of <T = (0.92±0.13)xl0-3mb. Although not clearly 
stated in their paper, this value does not correspond to 
the total cross section but represents the cross section for 
the transition to the first excited ^^F state only [l^l- As 
will be shown later, the data point of Becker et al. 3 
agrees with the data of Chow, Griffiths and Hall @ and 
Morlock et al. [7|. 



Hester, Pixley and Lamb 3 measured the total 
"'^^0(p,7)^^F cross section at six energies in the bom- 
barding energy range of i?p=140-170 keV. They mea- 
sured the ^^F activity after exposing an infinitely thick 
oxygen target to an intense proton beam. Thus they 
could not distinguish between transitions to the ground 
or first excited state in -'^^F. Their study represents an ab- 
solute cross section measurement, in the sense that the 
deduced cross section depends on absolute detection ef- 
ficiencies, absolute incident charge integration, stopping 
powers and target stoichiometry. Recall that many of the 
systematic errors related to these sources are avoided in 
the studies of Refs. d, 0] mentioned above. 

The cross section of Ref. 3 was obtained from the 
expression 



2£;3/2 



E 

z7 



(2) 



E. Other data 

Two more studies 0, Q report low-energy cross sec- 
tions for the ^^0(p,7)^''F reaction. These report total 
cross sections and provide no information on individual 
transitions to the two final "'^'^F states. We disregarded 
these data for the reasons given below. 

The study by Tanner @ reports cross sections in the 
bombarding energy range of Ep=274-616 keV. The ab- 
solute cross section measurement was performed with a 
WO3 target by using the activation method. At the time, 
stopping powers for tungsten were not available and, as 
the closest approximation, tabulated values for tantalum 
were used instead (23[. Clearly, these stopping powers 



are several decades old and the cross sections reported 
by Tanner Q should in principle be corrected by using 
modern stopping power values for oxygen and tungsten. 
However, our attempt at a correction was futile, mainly 
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because Tanner obtained the cross section from the yield 
and the effective stopping powers of his 150 keV thick tar- 
get by numerical integration. Since the measured yields 
are not provided in Ref. [3] there is no obvious way to 
correct his reported cross sections. Consequently, we are 
compelled to disregard the data of Tanner 5] in our anal- 
ysis. 

Rolfs [1| reports total ^^0(p,7)^^F cross sections in the 
bombarding energy range of Ep=0. 3-3.0 MeV. He did not 
determine the absolute scale of the cross section, but nor- 
malized his results relative to the absolute cross section 
of Tanner Q that was measured at an energy of Ej,=616 
keV. Since we did not succeed in correcting the latter 
data for improved stopping powers (see aboveV, we are 
also compelled to disregard the data of Rolfs 0] in our 
analysis. Nevertheless, the data of Ref. 0, which were 
measured over a broad energy range, demonstrated the 
rise of the S-factor at low energies and the measured 
energy dependence agrees with later results [Tj] . 



III. NUCLEAR REACTION MODEL 
CALCULATIONS 
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In this section we describe the fitting of a number of 
models to the data reported in Chow, Griffiths and Hall 
and Morlock et al. 7]. The fits are performed for 
each transition separately. No attempt has been made 
in the present work to constrain the fitting by using the 
reported total S-factors of Hester, Pixley and Lamb 
However, the latter data, together with the single data 
point of Becker et al. [1] for the transition to the first 
excited ^^F state, provide important cross checks for the 
present results. One of our main goals is to extrapo- 
late the S-factor to center-of-mass energies of E«50-140 
keV, corresponding to stellar temperatures of T«0.06- 
0.1 GK that are characteristic of hot bottom burning 
in AGB stars (see the shaded area in Fig. [2]). It will 
be seen that our procedures give rise to significantly im- 
proved thermonuclear reaction rates over the entire tem- 
perature range of astrophysical interest (T=0. 01-10 GK). 
Throughout this theoretical section and the following sec- 
tions, all kinematic quantities are given in the center-of- 
mass system, unless mentioned otherwise. 

A proper analysis should distinguish between statisti- 
cal and systematic errors. We corrected the data of Refs. 
0)0] for a number of different effects (Sec. |Tl|. However, 
the original papers do not provide sufficient information 
to quantitatively disentangle statistical and systematic 
errors. In view of this ambiguity, we decided to treat the 
combined uncertainties as statistical errors in the fitting 
procedure. The same problem and strategy is encoun- 
tered in Ref. [23 | . More information regarding our error 
analysis, and a justification for the adopted procedure, is 
given in Sec. HVl 

The radiative capture cross section for ^^0(p,7)^^F 



FIG. 2; (Color online) Total astrophysical S-factor for the 
^^0(p,7)^^F reaction versus center-of-mass bombarding en- 
ergy; (top) all data considered, (bottom) expanded view for 
low energies. A number of different corrections are applied to 
the data of Refs. 0, d, 0] while the results of Refs. [3, |^ are 
disregarded in the present work; see text for details. Note that 
we neglected the lowest energy data point of Ref. Q because 
the associated error is very large (±52%). The total S-factor 
varies slowly, except near the vicinity of the lowest lying res- 
onances at Ep=2.66 and 3.47 MeV The shaded area in 
the lower part indicates the effective energy range of stellar 
burning in massive AGB stars corresponding to temperatures 
of T=60-100 MK. 



from an initial state J^TTi to a final state J/tt/ is given by 



ae(J.x. ^ JfTTf) = — (2J, + l)|C/^(J,7r, ^ J/^/)P,(3) 



where U'^ is the transition matrix element. For a given 
final state Jy, the total cross section is obtained by sum- 
ming over all initial angular momenta J^. This definition 
of the cross section is common to all models; the only 
difference is the way of calculating the transition matrix 
element. In the potential model, this term is computed 
from an integral involving the initial and final wave func- 
tions, and the radial part of the electromagnetic operator. 
In the _R-matrix approach, it is split in two parts, involv- 
ing different parameters fitted to experiment. We give 
more detail below on both models. 
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A. Potential model 

We start the analysis of existing data by using a simple 
model for the low-energy cross section in ^^0(p,7)^''F. 
The model is referred to as "direct proton capture" and it 
assumes a single-step process where the proton is directly 
captured, without formation of a compound nucleus, into 
a final bound state with the emission of a photon. Many 
nonresonant reactions in nuclear astrophysics have been 
described by this model (see Ref. [1^ and Tab. I in Ref. 

M)- _ 

Numerous studies (see, for example, Ref. [26|) have 
shown that at relatively low bombarding energies and 
small binding energies of the final states the direct cap- 
ture of protons occurs mainly far outside the nuclear ra- 
dius. For this reason, the direct capture process is some- 
times referred to as "extra-nuclear" capture. This also 
implies that the calculated cross section is expected to 
be relatively insensitive to the details of the model used 
to describe the nuclear interior. 

In the potential-model formalism the transition matrix 
element is given by 



u,{r,E)0^xir)uf{r)rUr (4) 



where C^S denotes the spectroscopic factor of the final 
state [27[, and Oui\ is the radial part of the multipole 



operator for electromagnetic radiation of character lu (ei- 
ther electric or magnetic) and multipolarity A; Ui denotes 
the radial part of the partial waves of the initial scatter- 
ing state and u/ is the radial wave function of final bound 
state. The model adopted here is similar to the one em- 
ployed in Ref. Q. The full formalism is given in Refs. 

For example, for the ^^0(p,7)^^F reaction the tran- 
sition to the ground state proceeds predominantly via 
El radiation and angular momenta of £i — 1,3 (J^ = 
3/2-, 5/2-, 7/2") ^ £j = 2 [Jf = 5/2+) (DC-^0), 
while the transition to the first excited state at Ea;=495 
keV proceeds via El radiation and angular momenta of 
i^ = 1 iJ^ ^ l/2-,3/2-) if ^ (Jf = 1/2+) 
(DC^495). Other transitions, that is, those of Ml or 
E2 character, are negligible for the direct proton capture 
on 1^0. See also Ref. [s^. 

The scattering and bound state wave functions are gen- 
erated by an optical potential 

V{r)^^VF{r)+Uo-^L-S + Vc{r) (5) 
r ar 

where Vf denotes a folding potential and ^ is the folding 
potential strength parameter. Folding potentials have 
the major advantage that the geometry of the potential 
is fixed by the folding calculation. In other words, the po- 
tential can only be changed by a variation of the strength 
parameter ^. The (weak) spin-orbit potential term has 



the usual Thomas form and is characterized by the spin- 
orbit strength parameter ^^.o., while the Coulomb poten- 
tial Vc is given by a uniformly charged sphere of radius 
Rc- The folding potential is determined using the effec- 
tive nucleon-nucleon interaction adopted from the well- 
established DDM3Y parametrization. For the volume 
integral per interacting nucleon pair and the root-mean- 
square radius we find values of Ji^=525.93 MeV fm^ and 
''_F,rms=3.311 fm, respectively, while we also adopt the 
latter value for the Coulomb radius. For details regard- 
ing the folding procedure, see Refs. 0, [Ml- 

For generating the bound state wave function, the pa- 
rameter ^ is adjusted to reproduce the binding energies 
of the ground and first excited state, while the param- 
eter ^s.o. can be constrained by the energy splitting of 
the lowest 5/2+ and 3/2+ states in ^'^F. We find val- 
ues of ^==1.0976 and ^s.o.=— 0.1757 fm^. For generating 
the scattering state wave function, it was shown in Ref. 
[jj that good agreement between the experimental cap- 
ture data and the calculation is found by using a value 
of ^ «1.0 and the same value for S^s.o. as for the bound 
state calculation. This potential also describes reliably 
the ^^0(p,p)^^0 elastic scattering data at low energies. 

Once the folding potential strength parameter ^ for 
generating the scattering wave function has been fixed, 
the spectroscopic factor C^S is the only adjustable pa- 
rameter in the above model. Our strategy will be as 
follows. Initially, the potential parameter is held con- 
stant at 'C=l-0 and the data of Chow, Griffiths and Hall 
and Moriock et al. are fitted independently. The 
least-square fits provide for each transition and data set 
the corresponding value of C^S with an associated error. 
Afterward, the above procedure is repeated by system- 
atically varying the value of ^ in order to investigate the 
sensitivity of the extrapolated S-factor. 

Numerous studies have shown that the spectroscopic 
factors for the transitions to the ground and first excited 
states in ^^0(p,7)"'^^F are close to unity (see, for example, 
Ref. [1^). It must be emphasized, however, that we are 
here not at all concerned with the extraction of accurate 
spectroscopic factors from the measured low-energy cross 
sections. The magnitude of C^S will depend on the po- 
tential parameter ^, as is apparent from Eqs. [3] and [51 
In the present analysis, the spectroscopic factors are sim- 
ply regarded as intermediate results, or scaling factors, 
whose derived uncertainty determines partially the error 
in the extrapolated astrophysical S-factor. 

Our potential model does not account for the pres- 
ence of resonances. Therefore, we only consider the 
center-of-mass energy range below E=2.4 MeV for the 
least-squares fits described in this section. The fits 
to the experimental S-factor data for the transition to 
the ^^F ground state are shown in Fig. [31 The re- 
sults are C^S^DC 0) = 1.09±4.0% {xled='^-^) for the 
data of Chow, Griffiths and Hall ^ and C^S{DC 
0)=1.22±0.9% {xl^^=l.l) for the data of Moriock et al. 
[31 . The absolute normalization of these two data sets 
differs thus by 12%. Figure [H shows the fits to the ex- 
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TABLE I: Astrophysical S-factors for ^'^0(p,7)^^F at ener- 
gies of E — 90, 500, and 1000 keV. The results shown here are 
obtained by fitting the data of Ref. [7]] for the dominant tran- 
sition to the first excited state only. The S-factor obtained 
in the folding potential calculation of Ref. [3| is used as a 
reference, Siet- The parameter ^ refers to the strength of the 
scattering potential. 



E (kcV) 


90 500 1000 




Sr,{{E) (keVb) 


7.03 4.20 3.39 
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perimental S-factor data of the transition to the first 
excited state in ^^F. In this case, we find C^S{DC — + 
495)=1.05±1.7% (x^ed=l-5) for the data of Ref. [6] and 
C^SiDC 495)==1.02±0.5% (x^^^^^O.S) for the data of 
Ref. Q. The difference in the absolute normalization of 
the two data sets amounts to 3%. Note that the upturn 
in the S-factor at low energies is well understood. See, 
for example, Refs. @, 

A different choice for the parameter ^ of the scattering 
potential changes not only the magnitude, but also the 
energy-dependence, of the calculated direct capture S- 
factor. The sensitivity of the extrapolated S-factor to this 
parameter will be explored at a center-of-mass energy of 
E=90 keV which is located at the center of the energy re- 
gion important for the hot bottom burning in AGB stars 
(shaded area in Figs. [3] and ID). We quote in the follow- 
ing the results obtained for the dominant transition to 
the first excited state in ^''F (DC^ 495). For this sen- 
sitivity study we use the corrected experimental data of 
Ref. [7| with their statistical errors only (that is, exclud- 
ing the additional 10% overall uncertainty; see Sec. II. B). 
Similar results are obtained for the data of Ref. Q . For 
the reference calculation (using ^=1) we find a spectro- 
scopic factor of C^S = 1.02 and an S-factor of 5'rcf(90 
keV)= 7.03 keVb for the DC-^495 transition. For exam- 
ple, a variation in the potential strength parameter ^ by 
10% changes both the spectroscopic factor obtained from 
a least-squares fit of the data of Morlock et al. and the 
extrapolated S-factor at E=90 keV by about 5%. The re- 
sults of our sensitivity study for center-of-mass energies 
of E=^Q, 500 and 1000 keV are given in Tab. [B We also 
find that calculations with Woods-Saxon scattering po- 
tentials give similar results, even when the Woods-Saxon 
parameters are varied over a relatively broad range. 

It is interesting to note that a variation in the poten- 
tial strength parameter ^ influences the calculated energy 
dependence of the S-factor which, in turn, gives rise to a 
noticeable variation in the Xred value for the adjustment 
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FIG. 3; (Color online) S-factor for the proton capture into 
the ^^F ground state. The black and red lines represent fits to 
the data of Morlock et al. 0] and Chow, Griffiths and Hall 
respectively. The fines are obtained with a fofding potentiaf 
modef (Sec. III. A). The 68%-confidence bands are rather 
narrow and are not dispfayed for reasons of cfarity. Onfy data 
befow a center-of-mass energy of E=2.4 MeV are incfuded 
in the fits. The shaded area in the fower part indicates the 
effective energy range for the hot bottom burning in AGB 
stars. 



of the spectroscopic factor C'^S. A prerequisite for such 
a study is a large number of data points as provided by 
Ref. Q. For the folding potential we find a broad 
minimum near ^ « 0.9. The value of Xrcd increases by 
1 if ^ is changed by about 10%. Thus, the energy de- 
pendence of the S-factor confines the potential strength 
parameter to the range explored in Tab. IH 

In summary, a variation of the potential parameter ^ 
within a reasonable physical range (10%) changes the ex- 
trapolated S-factor by wl-5%, depending on the energy 
E. This value can be compared to the statistical error 
obtained from the least-squares fits alone (0.5-4%) and 
to the difference between the absolute normalizations of 
the two data sets under consideration (3-12%). 

The sum of the S-factors for both transitions is shown 
in Fig. [5] (top). The red and black lines indicate the fits 
to the data of Chow, Griffiths and Hall Q and Morlock 
et al. 0, respectively. It can be seen that the analyses of 
both data sets give consistent results. As a cross-check, 
the low-energy data of Hester, Pixley and Lamb [3] are 
also shown. Recall that these authors only report the 
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FIG. 4; (Color online) S-factor for the proton capture into 
the first excited state of ^^F. The black and red lines represent 
fits to the data of Morlock et al. ^Ti] and Chow, Griffiths and 
Hall 01, respectively. The lines are obtained with a folding 
potential model (Sec. III. A). Only data below a center-of- 
mass energy of E=2.4 MeV are included in the fits. Notice 
the single data point reported by Becker et al. shown in 
green. See also caption of Fig. [S] 
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FIG. 5: (Color online) Total S-factor for ^^0{p,y)"F at low 
bombarding energies. The red and black lines indicate fits to 
the data of Chow, Griffiths and Hall @| and Morlock et al. 
3, respectively, by using a folding potential model (top part; 
see Sec. HI. A) and an R-matrix model (bottom part; see Sec. 
III.B). The low-energy data of Hester, Pixley and Lamb [3] 
are shown for comparison. 



total S-factor. Clearly, the solid lines are consistent with 
the data of Ref. 4], providing further support to our 
results. 



B. R-matrix model 

1. Formalism 

The R-matrix theory combines the flexibility of a phe- 
nomenological approach with the physical content of the 
Coulomb wave functions. It can be applied to resonant as 
well as to non-resonant reactions. In the R-matrix theory 
the relative coordinate r between the colliding nuclei is 
divided into two sectors, limited by the channel radius a. 
In the internal region (r < a) the physics of the problem 
is determined from the properties of poles. In the exter- 
nal region (r > a) , the relative wave function is given by 
a linear combination of Coulomb wave functions. This 
procedure provides an accurate extrapolation down to 
low energies, where Coulomb effects are expected to play 
a major role. The poles are associated with resonances or 
bound states of the unified nucleus, and their properties 



(energy, particle and 7-ray partial widths) are fitted to 
the available data. 

The formalism of the R-matrix method is well known 
[S^l- Here we use the formulation of Ref. The 
electromagnetic transition matrix element If^ involved 
in Eq. ([3|) is divided into two parts, associated with the 
internal and external regions, as 

[/^ = Ul, + Ul, (6) 

We present the different contributions in Ref. [24]. The 
internal part U^^^ is given by the energies and partial 
widths of the poles, while the external contribution C/J^j 
is provided by an integral involving Coulomb wave func- 
tions and is evaluated from the channel radius a to infin- 
ity. This external term is proportional to the asymptotic 
normalization constant (ANC) of the final (bound) state. 
Let us point out that both contributions are related to 
each other, as the external term involves the phase shift 
which, in turn, depends on the R matrix. 

The R-matrix theory can be easily extended to non- 
resonant reactions. In that case one usually assumes that 
the non-resonant contribution is simulated by a high- 
energy background pole. Although this method is accept- 
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able, it presents some drawbacks: the results are some- 
what sensitive to the choice of the pole energy, which in 
addition has no genuine physical meaning. This prob- 
lem can be addressed by assuming that the R-matrix is 
constant. Non-resonant transfer reactions have been in- 
vestigated previously in that way [33|. The R-matrix is 
therefore given by 

R{E) = Rl (7) 

where subscript p refers to the proton channel. When 
-Rq = 0, the problem is equivalent to the hard-sphere 
approximation. 

This procedure can be extended to capture reactions. 
Let us consider the internal contribution in Eq. ([6]). Ac- 
cording to Ref. [ill, we have in a general case 



1 ^ 

U^nt = «^ exp(zfes) ^ _^^ ^ 

i—1 



where i?g is a constant connected with the electromag- 
netic transition. This parametrization is equivalent to 
the usual "background" term, but without the redun- 
dancy associated with the energy of the background pole. 
Note that a non-resonant capture transition is character- 
ized by the fitting parameters i?Q, i?Q and the ANC of 
the final state. 



2. Results 

For the ground-state transition in -'^^0(p,7)"'^^F, we 
have ^i = 1 [Ji = 3/2^), and a small contribution 
from 4 = 3 (J, = 5/2-, 7/2-). All these transi- 
tions are non-resonant. For the excited-state contribu- 
tion {Jf = 1/2+), we include the l/2~ resonance at 
Er — 2.50 MeV, with parameters taken from Ref. 
The Ji — 3/2^ partial wave is also included as a non- 
resonant contribution. For simplicity we assume that 
the partial waves of the same spin-orbit doublet (that 
is, Ji = 5/2~ and Ji = 7/2~) have the same R-matrix 
parameters. For the channel radius we used a value of 
a=5 fm. A simultaneous fit for the ground and excited 
states is then performed. 

The R-matrix fits for the transitions to the ^^F ground 
and first excited state are shown in Figs. [B]and[71 respec- 
tively. For the asymptotic normalization constant of the 
ground state transition we find values of ANC=1.19±0.02 
fm-i/2 [xled=^-'^) for the data of Chow, Griffiths and 
Hall and ANC=1.13±0.01 fm-^/^ {xl^^=O.Q) for the 



where TV is the number of poles, F^, F^^i) are their 
energies, particle and 7-ray partial widths, (5hs is the 
hard-sphere phase shift, and the constant L is related 
to the Coulomb wave functions (see Refs. [H, H?] for 
details). Including the energy dependence of the widths 
we have 



f J^,, = 2P,{E)^^f°^^,{E - Eff^+' (9) 



where Ef is the energy of the final state, A is the order of 
the multipole, and jf, r° ^ are energy-independent quan- 
tities. For a non-resonant reaction {N = 1), Eq. ^ can 
be rewritten as 



(10) 

I 

data of Morlock et al. [7|]. The absolute normaliza- 
tion of these two data sets differs by 5%. For the 
transition to the first excited ^^F state, the results are 
ANC=81.0±0.9 fm"i/2 for the data of Ref. 

and ANC=82.3±0.3 fm^^/^ {xl^^=0.7) for the data 
of Ref. 0, where the absolute normalization differs by 
2%. Note that the asymptotic normalization constants 
from our R-matrix analysis are consistent with the mi- 
croscopic results of Ref. [sj . 

The sum of the S-factors for both transitions is dis- 
played in Fig. O (bottom). The red and black lines show 
the results from fitting the data of Chow, Griffiths and 
Hall and Morlock et al. , respectively. It is apparent 
that the R-matrix analyses of both data sets give consis- 
tent results and describe the total S-factor data rather 
well. The channel radius is adopted here as a = 5 fm, 
a typical value used in the literature [23|. The sensitiv- 
ity of the extrapolated S-factor to systematic variations 
of the channel radius a were investigated at a center-of- 
mass energy of 90 keV. We find that a variation of the 
channel radius by 10% changes the extrapolated S-factor 
at E=90 keV by about 2%. 



IV. RECOMMENDED S-FACTOR AND 
REACTION RATES 

In the previous sections we presented the total S- 
factors and their components for two different data sets, 
analyzed by using two different nuclear reaction models. 



I 

Ul, = ^' exp(z<5Hs)T3^ [2<i?2A. iE){E Efr+'] ' 
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FIG. 6; (Color online) S-factor for the proton capture into 
the ^^F ground state. The black and red lines represent fits to 
the data of Morlock et al. 0] and Chow, Griffiths and Ifall ^] , 
respectively. The lines are obtained with an R-matrix model 
(Sec. III.B). The shaded area in the lower part indicates the 
effective energy range for the hot bottom burning in AGB 
stars. 



FIG. 7: (Color online) S-factor for the proton capture into the 
first excited state of ^^F. The black and red lines represent 
fits to the data of Morlock et al. |7| and Chow, Griffiths 
and Hall [^, respectively. The lines are obtained with an 
R-matrix model (Sec. III.B). The shaded area in the lower 
part indicates the effective energy range for the hot bottom 
burning in AGB stars. 



We will now describe our method of estimating recom- 
mended S-factors and their associated errors. This in- 
formation will then be used in order to calculate new 
recommended reaction rates. 

The method is explained in Fig. [5] showing, as an ex- 
ample, the total S-factor at an energy of E=0.090 MeV, 
that is, at the center of the effective energy range im- 
portant for the hot bottom burning in AGB stars. The 
left and middle panels present results obtained from the 
folding potential model and the R-matrix approach, re- 
spectively. The x-axis labels "CH" and "MO" refer to 
the analysis of the data of Ref . [Q] and Ref. [7] , respec- 
tively. The solid bars indicate la errors resulting from 
the fitting of the data alone (statistical errors). Note that 
in general the errors obtained in the R-matrix approach 
with three fitting parameters in our case are larger than 
those derived from the potential model with just one fit- 
ting parameter. The errors indicated by the open bars 
are obtained if systematic model uncertainties (i.e., those 
obtained by varying the parameters of the model within 
a reasonable physical range) are added quadratically to 
the statistical errors. The numerical values of the statis- 
tical and total errors (in percent) are listed at the top 
and bottom, respectively, of each error bar. The label 



"AV" refers, for each model, to the weighted average of 
the two data sets, where only statistical errors are used 
for the weighing procedure and systematic model errors 
are added quadratically afterwards. The right panel dis- 
plays the final adopted S-factor (label "AD"). Its total 
error is given by the extrema of the average values re- 
sulting from the two reaction models. The mean value 
of the final recommended S-factor is then given by the 
arithmetic average of the upper and lower bound of the 
error bar (open bar in right panel) . 

It was pointed out in Sec. Ill Bl that an additional un- 
certainty of 10% has to be taken into account for the 
data of Morlock et al. Unfortunately, insufficient infor- 
mation is provided in Ref. 0] to determine how much of 
this uncertainty arises from systematic and from statis- 
tical effects. Recall that in the above analysis we treat 
the entire additional uncertainty as a statistical error in 
the fitting procedure (see Sec. IIII|1 . Our treatment will 
tend to underestimate the true error. Alternatively, one 
may regard the entire additional uncertainty as a sys- 
tematic error, which would tend to overestimate the true 
error. Clearly, the true error must be located somewhere 
between the two extremes. Test calculations have been 
performed in order to quantify this effect. Consider as 
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FIG. 8: Total S-factor for ^"^0(p,7)^^F at a center-of-mass 
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present and the NACRE results, respectively. 



an example the folding potential analysis of the domi- 
nating transition to the first excited state for an energy 
of E=0.090 MeV (similar to the left-hand side of Fig. 
8) . The alternative assumption of treating the additional 
uncertainty as a systematic error would increase the av- 
eraged ("AV") S-factor by 2.6%, while its uncertainty 
would increase from 5.1% to 5.3%. Similar results are ob- 
tained in the R-matrix analysis. It is obvious that such 
small variations are negligible for the adopted S-factor 
("AD"). 

The thermonuclear rate per particle pair for a reaction 
involving two nuclei is given by 

NA{av) = N^^^^^ [ S{E) e-^/'^^-- dE 

(11) 

where S{E) is the S-factor at energy E and the factor 
^~E/kT (derives from the Maxwell-Boltzmann distribu- 
tion. The quantities T, k, Na and m are the plasma 
temperature, the Boltzmann constant, Avogadro's con- 
stant and the reduced mass of the interacting nuclei, re- 
spectively. The thermonuclear rates for the ^^0(p,7)^^F 
reaction, calculated by integrating Eq.[Tl] numerically us- 
ing the S-factor recommended in this work, are listed in 
Tab, [m for temperatures in the range of T=0.01-2.5 GK. 
For higher temperatures, the reaction rates are influenced 
by energies not covered in the present work (E>2.4 MeV). 

Results for the reaction rates of ^^0(p,7)^^F are shown 
in Fig. [9] The curves display the ratios of lower bound, 
recommended rate and upper bound on the rate over the 
present recommended rate. Solid and dashed lines in- 
dicate the rate ratios resulting from our work and from 
NACRE 0, respectively. It can be seen that the ratio of 
the NACRE and the present recommended rates is close 
to unity. However, below a temperature of T=0.5 GK, 
the present reaction rate errors amount to only ^7%. 



This represents an improvement by a factor of w4 over 
the results published by NACRE which followed a rather 
conservative rate evaluation procedure. Note that be- 
cause of extensive efforts, the ^^0(p,7)^^F reaction has 
now the most precisely known rate involving any target 
nucleus in the mass A > 12 range. Astrophysical conse- 
quences for the oxygen isotopic ratios from hot bottom 
burning in AGB stars are presented in the following sec- 
tion. 



V. HOT BOTTOM BURNING IN AGB STARS 

The aim of this section is to test the effect of our new 
^^0(p,7)^^F reaction rate on the ^^O/^^O ratios pro- 
duced by HBB in massive AGB stars. This is of rele- 
vance because of the presence of a very small fraction of 
stellar meteoritic oxide grains that might have originated 
in massive AGB stars. These stars may have contributed 
to the abundances of some radioactive nuclei in the early 
solar system, for example, ^^Al and ^"Fe, as well as to 
other abundance anomalies observed in meteoritic mate- 
rial, for example, ^^Ni, ®^Rb, and ^^Zr [3^. It is impor- 
tant to verify if there is also a population of meteoritic 
stellar grains that may have originated from this stellar 
site. 

The few (in fact, four) grain candidates for a massive 
AGB origin exhibit number abundance ratios of ^^O/^^O 
< 10-^, i^O/i^O between 7 x 10^* and 1.5 x lO'^, and 
large ^^Mg excesses caused by the decay of ^^Al, with 
inferred initial ^SAl/^^Al ratios of ~ 0.01-0.09. One of 
them, spinel grain 0C2, also shows an excess in ^^Mg. 
Such signatures could be produced in massive AGB stars 
via the combined activation of HBB and the operation of 
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the ^^Ne+a reaction during unstable He burning (ther- 
mal pulses). For a detailed discussion, see Lugaro et al. 
[T^ l. We focus here on the ^^O/^^O ratio, which is the 
most problematic signature to be interpreted as produced 
by HBB (see also Ref. [36.] )• 

During HBB the abundance of ^^O quickly reaches 
equilibrium. Thus the ^^O/^^O ratio is mostly deter- 
mined by the ratio of the reaction rates that produce and 
destroy "O, that is, the i60(p,7)i^F and i^O(p,a)i'*N 
reactions, respectively. Because the laboratory measure- 
ments of the composition of meteoritic stellar grains pro- 
vide isotopic ratios of very high precision, with error bars 
on the order of a few percent only, small uncertainties in 
the reaction rates are indispensable in this case to address 
the possible match of the models with the observations. 

We consider here a stellar model of a typical massive 
(6.5 Mq) AGB star of solar metalhcity {Z = 0.02) expe- 
riencing HBB. This model was computed by Karakas & 
Lattanzio ^32] using mass loss on the AGB from Ref. (sH] 
and it experienced 39 thermal pulses. The temperature 
at the base of the convective envelope during the AGB 
phase increases from T=60 MK to T=87 MK at the 25**^ 
thermal pulse, and then decreases again to T=20 MK at 
the end of the computed evolution. The high temper- 
ature reached by this model allows for the production 
of enough ^^Al to match the composition of grain 0C2, 
while the ^^O/^^O ratio could be reproduced within the 
previous M "'^^0(p,7)^^F reaction rate uncertainties (see 
Fig. 7 of 111). 

Similar to the procedure of Lugaro et al. [T3l |. we an- 
alyze the impact of the reaction rate uncertainties using 
a post-processing nucleosynthesis code. Hence we disre- 
gard any feedback of the different reaction rates on the 
stellar structure. This is justified in this case because: 
(i) the ^^0(p,7)^''F and ^^0(p,a)^^N reactions generate 
only a small amount of energy within the CNO cycle; 
and (ii) the ^^0(p,7)^^F reaction is only marginally ac- 
tivated in our model (only 23% of the initial ^^O is de- 
stroyed by the second dredge-up and by HBB). For the 
^^0(p,7)^^F reaction we use our new rates, while for the 
^^0(p,a)^'*N reaction we employ the same rate as in Ref. 
[Tl which is based on the latest available experimental 
information [H, H^. Different runs of the 6.5 Mq stel- 
lar model are then computed by varying independently 
the rates of these two reactions between their lower and 
upper bounds. The results of this procedure are given in 
Tab. IHIl showing how the errors in the rates propagate 
to uncertainties of the computed '^^O/^^O ratio. Note 
that the calculated ^^O/^^O ratio is almost the same as 
in Ref. [ll because the recommended ^^0(p,7)^^F reac- 
tion rate did not change significantly (see Fig. [9]) . The 
main result here is that the total uncertainty in the de- 
rived abundance ratio is reduced from -1-63% and —73% 
when using the NACRE rate for if^0(p,7)"F to -f35% 



and —30% when employing our new reaction rate. 

In conclusion, the measured -'^^O/^^O ratio of grain 
0C2 (= 1.25 ± 0.07 X 10-3) could be reproduced 
within the large error bars of the NACRE compilation 
(2.44+i[ 7g X lO-'^) in models of massive AGB stars; how- 
ever, the much more precise ^^0(p,7)^^F rate of the 
present work leads to 2.52^1^:^^ x lO'^ for the "O/i^Q ra- 
tio and disagrees with the measured value. Consequently, 
there is not clear evidence to date for any stellar grain 
origin from massive AGB stars. Stellar model uncertain- 
ties, such as different mixing prescriptions ^4.1] and mass 
loss rates [i^ still need to be carefully evaluated in this 
context. Another possibility is that we have not yet dis- 
covered grains from massive AGB stars because they are 
perhaps smaller in size than t he g rains currently analyzed 
in the laboratory (> 1 /xm) |43l 14^. Clearly, a solution 
to this problem requires future efforts. 



VI. SUMMARY 

In this work we focused our attention on a reanaly- 
sis of the ^^0(p,7)^^F reaction rates. We started from 
an evaluation of all the original data and performed a 
number of corrections when appropriate. The modi- 
fied and improved data are then interpreted in terms of 
two independent models of nuclear reactions, a poten- 
tial model and an R-matrix approach. This attention 
to detail was clearly not practical in previous reaction 
rate evaluations. After combining the results from the 
two reaction models, we find ^^0(p,7)^^F reaction rate 
errors of less than 7% over the entire range of astrophys- 
ical interest. In other words, the ^^0(p,7)^^F reaction 
exhibits now the most precisely determined thermonu- 
clear rates among any charged-particle capture reactions 
in the A > 12 mass range. Our new results are then 
incorporated into models of massive AGB stars in order 
to study the derived oxygen isotopic ratios. Contrary to 
previous conclusions, we find now that there is presently 
no clear evidence of a massive AGB star origin for any 
observed stellar grains. 
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TABLE II: Total thermonuclear reaction rates for 
^''0(p,7)"F in units of (cm^mor^s-^). 



T (GK) 


Low 




Recommended 


High 


0.01 


6.674x10" 




7.200x10"^^^ 


7.733x10"^^ 


0.011 


7.026x10- 




7.578x10"^* 


8.138x10"^* 


0.012 


5.638x10- 


— 23 


6.080x10"^'^ 


6.528x10"^-^ 


0.013 


3.626x10- 




3.910x10"^^ 


4.197x10"^^ 


0.014 


1.941x10 




2.093x10"^! 


2.246x10"^! 


0.015 


8.910x10- 




9.604x10"^! 


1.031x10"^° 


0.016 


3.587x10- 




3.866x10"^° 


4.148x10"^° 


0.018 


4.211x10- 




4.536x10"!'' 


4.866x10"!® 


0.02 


3.505x10- 


— 18 


3.775x10"!^ 


4.048x10"!® 


0.025 


2.431x10- 




2.616x10"!*^ 


2.803x10"!'' 


0.03 


6.124x10- 




6.586x10"!^ 


7.053x10"!^^ 


0.04 


6.659x10- 


— 13 


7.155x10"!^ 


7.657x10"!^ 


0.05 


1.847x10- 




1.984x10"!! 


2.123x10"!! 


0.06 


2.309x10- 




2.481x10"!° 


2.655x10"!° 


0.07 


1.726x10- 




1.855x10"°^ 


1.985x10"°® 


0.08 


9.029x10- 


—09 


9.706x10"°^ 


1.039x10"°* 


0.09 


3.644x10 




3.917x10"°* 


4.193x10"°* 


0.1 


1.208x10- 




1.299x10"°'^ 


1.391x10"°^ 


0.11 


3.439x10 


—07 


3.697x10"°^ 


3.957x10"°'' 


0.12 


8.662x10- 


—07 


9.312x10"°'' 


9.966x10"°'' 


0.13 


1.975x10- 




2.124x10"°*^ 


2.273x10"°® 


0.14 


4.149x10- 




4.460x10"°^ 


4.773x10"°® 


0.15 


8.133x10- 


—06 


8.742x10"°^ 


9.355x10"°® 


0.16 


1.504x10- 




1.616x10"°^ 


1.729x10"°^ 


0.18 


4.449x10- 




4.781x10"°^ 


5.114x10"°^ 


0.2 


1.129x10- 


—04 


1.213x10"°* 


1.298x10"°'! 


0.25 


7.231x10 


—04 


7.766x10"°* 


8.305x10"°* 


0.3 


2.954x10- 


-03 


3.172x10"°^ 


3.391x10"°^ 


0.35 


9.035x10- 


-03 


9.696x10"°^ 


1.036x10"°^ 


0.4 


2.262x10- 


-02 


2.426x10"°^ 


2.592x10"°^ 


0.45 


4.896x10- 


-02 


5.250x10"°^ 


5.605x10"°^ 


0.5 


9.495x10- 


-02 


1.018x10"°! 


1.086x10"°! 


0.6 


2.813x10 


-01 


3.012x10"°! 


3.212x10"°! 


0.7 


6.659x10 


-01 


7.122x10"°! 


7.587x10"°! 


0.8 


1.350x10+°° 


1.442x10+°° 


1.535x10+°° 


0.9 


2.447x10+°° 


Z.OiOX iU 


o 'y'yc^ 1 n"t~00 


1.0 


4.074x10+°° 


4.340x10+°° 


4.607x10+°° 


1.25 


1.123x10+°^ 


1.192x10+°! 


1.261x10+°! 


1.5 


2.418x10+°! 


2.557x10+°! 


2.696x10+°! 


1.75 


4.437x10+°! 


4.677x10+°! 


4.918x10+°! 


2.0 


7.295x10- 


fOl 


7.667x10+°! 


8.041x10+°! 


2.5 


1.586x10+°^ 


1.659x10+°^ 


1.731x10+°^ 
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TABLE III: Final ^^0/^*^0 abundance ratios for 6.5Mo, Z = 
0.02 stellar models computed with all possible combinations 
of ^''0(p,7)^''F and "0(p,q)^''N reaction rates: LR^lower 
bound of rate (at T = 60 — 100 MK), RR=recommended rate, 
and UR=upper bound on rate (at T = 60 - 100 MK). The 
17q^16q j-g^^jQ jg explicitly given for the RR-RR combination 
in the center of the table, while percent variations with respect 
to this case are given for the other rate combinations. Results 
for the ^^0(p,7)^'^F rate derived in this paper are in roman 
font, results using the NACRE % rate are in italics. 





LR (-7%, 


-13%) 


^«0(p,7)^'F 
RR 


UR (+7%, +30%) 


LR (~ -22%) 


+17%, - 




+26%, +25% 


+35%, +63% 


i^O(p,a)"N RR 


-7%, - 


43% 


2.52xl0"^ 2.44x10-'^ 


+7%, +30% 


UR (~ +18%) 


-30%, - 


-75% 


-20%, -21% 


-13%, +7% 



